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Motivation 


Migmtiom The process by which a propagating delamination 
relocates to a new ply interface via matrix cracking 


Impact 



M. McElroy et al. A numerical and 
experimental study of damage growth in 
a composite laminate , in proceedings of 
the ASC 29th Technical Conference, 
San Diego, CA, USA, 2014. 


Skin-stringer pull off 


R. Krueger et al. Fatigue Life 
Methodology for Bonded Composite 
Skin/Stringer Configurations. 
NASA/TM-2001 -21 0842, 2001 . 
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Experiments: delamination migration 
Test Setup - Premise 


Delamination 
(“positive” shear stress) 


Migration 

(“negative” shear stress) 



*E.S. Greenhalgh, C. Rogers, P. Robinson. Fractoaraohic observations on delamination growth and the subsequent migration 
through the laminate . Composites Science and Technology, 69:2345-2351 , 2009. 






J.G. Ratcliffe, M.W. Czabaj and T.K. O’Brien. A test 
013 - 218028 , 2014 . 


haracterizing delamination migration in carbon/epoxv tape laminates 


Experiments: delamination migration test 
Test setup 

Cross-ply laminate 
“2D” migration process 
Pre-crack (Teflon insert) 
between 0° and 90° ply 12.7 a °" 49 

Variable load position (L) . 
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Experiments: delamination mi 
Test setup - overview 


Delamination 
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Experiments: delamination migration test 
Test setup - overview 







Experiments: delamination migration test 
Test setup - validation data 



Damage morphology 



Load - displacement 



Migration location 
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Experiments: delamination migration test 


© Modeling approach: Floating Node Method (FNM) and 
Virtual Crack Closure Technique (VCCT) 
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Floating Node Method 


Same implementation 
strategy suitable for standard 
finite element architecture 
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Method (PNM) 
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Floating Node 
Method (FNM) 
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Remeshing 
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Same solution 


Same solution 


B.Y. Chen, S.T. Pinho, N.V. De Carvalho, P.M. Baiz, T.E. Tay, A floating node method for the modelling of discontinuities in 
composites , Engineering Fracture Mechanics, Vol. 127:104-134, 2014. 


Floating Node Method (FNM) 
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Floating Node Method (FNM) 
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Floating Node Method (FNM) 


T crack Intersecting cracks 



Key Characteristics: 

• Floating Nodes are topologically related to each element with no initial position 
assigned 

• The position of the floating nodes is assigned only after the crack path is determined 

• The floating nodes are used to form sub-elements within the original element and 
accommodate crack networks 

• Ideally suited to represent multiple cracks and their intersection 

• Can be coupled with Virtual Crack Closure Technique (VCCT) and cohesive zone 
crack formulations to model crack propagation 
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Floating Node Method & Virtual Crack Closure Technique 


Virtual Crack Closure Technique (VCCT): 



Mode I 
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Mode II 
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FNM & VCCT applied to cross-ply laminates: 


Laminate 

[ 0 °/ 90 ° 2 / 0 °] 

f 1 0° 


90 ° 


0 ° 


N.V. De Carvalho et al, Modeling delamination migration in cross-ply 
tape laminates. Composites Part A: Applied Science and 
Manufacturing, 71, 192-203, 2015. 


1 FNM Element 

(multiple plies) 



A Floating node (DoF) 

□ Coordinates of crack 
positions 
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FNM & VCCT applied to cross-ply laminates: 


Laminate 

[ 0 °/ 90 ° 2 / 0 °] 

f 1 0° 


90 ° 


0 ° 


1 FNM Element 

[ 0 °/ 90 ° 2 / 0 °] 



• Real node 
A Floating node (DoF) 

□ Coordinates of crack 


positions 
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FNM & VCCT applied to cross-ply laminates: 


Quasi-static 

• Fracture Criterion: 

/ {GuG n ) = -1 = 0 

• Mixed Mode exponential law: 

Gi nt = Gj c + (Guc — Gj c ) ^ 

Fatigue 

— A ( Grmax) n 

. / \ ( Gumax \ 

n = ri! + (jin ~ ni) I — J 

A = Aj + (A n - A:) 


Delamination 



# Real node 
A Floating node (DoF) 

□ Coordinates of crack 
positions 
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FNM & VCCT applied to cross-ply laminates: 
Migration onset 


Quasi-static 


winter 

Lt c 


> l 



Fatigue 
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VdJV ), F t > 0 


Migration onset 

(delamination to matrix crack) 



# Real node 
A Floating node (DoF) 

□ Coordinates of crack 
positions 
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FNM & VCCT applied to cross-ply laminates: 
Migration onset - quasi-static 
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FNM & VCCT applied to cross-ply laminates: 
Migration onset - quasi-static 
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FNM & VCCT - application to c< 
Migration onset - fatigue 
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FNM & VCCT applied to cross-ply laminates: 


Quasi-static 


Matrix Crack 


/ (Gj, Gu) = —f — 1 = 0 

Cr Ic 

Fatigue 

— — : = Aj ( GTmax) ni 
diV 


Maximum tangential stress criterion: 




# Real node 
A Floating node (DoF) 

□ Coordinates of crack 


positions 
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FNM & VCCT - application to composites: migration matrix 
crack to delamination interaction 


• Topological criterion 
- local delamination is 
onset when matrix crack 
reaches interface 


Migration 

(matrix crack to delamination) 



# Real node 


A Floating node (DoF) 

□ Coordinates of crack 
positions 
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Fatigue algorithm 



DETERMINE THE GROWTH 
W RATE FOR EACH CRACK 


DETERMINE THE NUMBER OF 
CYCLES NEEDED TO 
Q PROPAGATE EACH CRACK BY 
ONE ELEMENT, AND THE 
CRACK WHICH PROPAGATES 
IN FEWEST CYCLES 

O PROPAGATE THE CRACK 


O ACCUMULATE THE CYCLES 








Verification - Static: DCB 


Load, 

N 



* R. Krueger. An Approach to Assess Delamination Propagation Simulation Capabilities in Commercial Finite Element Codes . NASA/TM-2008- 
215123, 2008 


Verification - Fatigue: DCB benchmark 



30 — 

1 100 10000 1000000 

N, Cycles 


* R. Krueger. Development of a Benchmark Example for Delamination Fatigue Growth Prediction . NASA/CR-201 0-21 6723 
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Experiments: delamination migration test 

Modeling approach: Floating Node Method (FNM) and 
Virtual Crack Closure Technique (VCCT) 


(3) Validation: modeling delamination migration 
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Validation: Delamination migration test 
Numerical model 



1 Stacking sequence: 

[904/03/(90/0) 2s /02/0/90 4 /T/0/904/0 2 /(90/0) 2s /0 2 /90 3 /0/90] 


I I 


II 


Plies modeled within the FNM 
element 


Model details 


Contact modeled between specimen and 
clamps/baseplate 

Clamping force applied in a first static step 
Abaqus/Standard (Implicit) + UEL 
All material properties obtained using 
standard/recommended test methods 
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1 FNM 
Element 


Dimensions (mm) 

~B* 2h C S ip 

12.7 5.25 12.7 115 49 


*B is the width of the specimen (out-of-the page); 

90° - specimen width direction; 0° - specimen span direction 
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Validation: delamination migration test 
Results - migration process 



Delamination 


Migration 


Observations 


• Correct sequence of events: delamination followed by migration 

• Failure morphology well captured - including crack path 
through-thickness 
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Validation: delamination migratio 
Results - load vs displacement 


L=1.0a,>: 



400 


300 


Load, 200 
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100 
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EXPERIMENTS 


— SIMULATIONS 


MIGRATION 


Migration 



0.5 1 .0 

Displacement (mm) 


Observations 


• Max load: good agreement 

• Delamination: unstable growth 
followed by arrest and 
subsequent unstable and stable 
growth 

• Migration: predicted before 
delamination arrest 
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Validation: delamination migration test 
Results - load vs displacement 


L=1.1a 



Load 
N 




Validation: delamination migration test 
Results - load vs displacement 


L=1.2a 



Load 

N 



Displacement, mm 


Validation: delamination migration test 
Results - load vs displacement 



Observations 


Load 
N 



Displacement, mm 
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Validation: delamination migration test 
Results - Migration location 
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Observations 


Trend well captured 
Conservative predictions 


1.2 


1.3 
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Fatigue - Preliminary results 
Delamination growth and cycles to migration 


Constant amplitude, R = 0.1 and f = 5 Hz: 



50,000 100,000 150,000 200,000 

N, Cycles 
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Summary 


• Developed a finite element model based on the 

Floating Node Method combined with the Virtual Crack 
Closure Technique to capture the interaction between 

delamination and matrix-cracking 

• Identified and applied migration criteria for both quasi- 
static and fatigue loading 

• Compared simulations and experiments. 

- Good agreement observed for load-displacement, migration 
location and path 

• Validation of the fatigue simulations are in progress 
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*R. Krueger. Development of and application of benchmark examples for mixed-mode I/ll quasistatic delamination propagation predictions. 
NASA-CR-2012-217562, 2012. 



